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Abstract 

We study the thermal stability of the periodic (P) and quasi-periodic (Q) precessional modes of 
the nanoparticle magnetic moment induced by a rotating magnetic field. An analytical method for 
determining the lifetime of the P mode in the case of high anisotropy barrier and small amplitudes 
of the rotating field is developed within the Fokker-Planck formalism. In general case, the thermal 
stability of both P and Q modes is investigated by numerical simulation of the stochastic Landau- 
Lifshitz equation. We show analytically and numerically that the lifetime is a nonmonotonic 
function of the rotating field frequency which, depending on the direction of field rotation, has 
either a pronounced maximum or a deep minimum near the Larmor frequency. 

PACS numbers: 75.50.Tt, 76.20. +q, 05.40.-a 



1 



I. INTRODUCTION 



Magnetic nanoparticles are of great interest because of their nanoscale physical properties 
and many current and potential applications. These applications range from high-density 
storage medial and spintronic devices^ to biomedical applications like drug delivery, cell 
separation, cancer treatment and many others (for a review, see Refs.[5|-[7(). Since the physical 
properties of nanoparticles play a decisive role in all these applications, their study is of 
fundamental importance. In particular, for high-density storage media, e.g., bit-patterned 
medial where each nanoparticle is a carrier of information, the thermal stability of a given 
direction (or magnitude) of the nanoparticle magnetic moment is one of the most important 
problems. The reason is that under thermal fluctuations the magnetic moment can be 
randomly switched to a new state leading to the loss of information. 

In the case of ferromagnetic nanoparticles, the fluctuation dynamics of the nanoparticle 
magnetic moment can be described by the stochastic Landau-Lifshitz equation. If the noise 
term in this equation is approximated by the Gaussian white noise, then the probability 
density of the magnetic moment satisfies the Fokker-Planck equation whose properties are 
well known.- This approach, introduced by Brown^ almost five decades ago, has become 



11 



and references 



an important tool in the study of stochastic magnetic dynamics (see Ref. 
therein). To characterize the thermal stability in the case of uniaxial nanoparticles, it is 
often enough to determine the lifetimes of the nanoparticle magnetic moment in the "up" 
and "down" states. In the above approach, the lifetime in a given state is usually interpreted 
as the relaxation time. However, from a theoretical point of view, the lifetime is reasonable 
to associate with the mean first-passage time (MFPT), i.e., average time that a random 
process dwells in a prescribed state. An additional advantageous feature of this definition of 
the lifetime is that the MFPT method is mathematically well developed.— - — This approach 
was first applied to study the magnetic relaxation in systems of noninteracting^ and dipolar 
interacting^ nanoparticles subjected to a constant magnetic field. 

In general, the lifetime depends on both intrinsic properties of nanoparticles and external 
magnetic fields. The case when the external fields contain a rotating magnetic field applied 
perpendicular to the easy axes of nanoparticles has a particular interest. On the one hand, 
this is because the rotating field plays a key role in the microwave-assisted switching^— and, 
on the other hand, because the corresponding dynamical equations (without accounting the 
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thermal fluctuations) can often be solved analytically.— — Specifically, it has been shown 
(see also Ref. 



281 ) that the rotating field can induce two types of the stable precessional modes 
of the magnetic moment, namely, the periodic (P) and quasi-periodic (Q) modes. Under 
certain conditions,— 1 ^ one mode can exist in the up state of the magnetic moment and the 
other in the down state. The thermal fluctuations make the random transitions between 
these modes possible, and the problem of the lifetime of the P and Q modes appears. 
Some aspects of this problem have already been considered previously in the context of 
magnetic relaxation and induced magnetization.—"— But the dependence of the lifetime 
on the parameters of the rotating field has not been studied systematically. At the same 
time, the effect of strong dependence of the lifetime on the rotating field frequency, which is 
expected to exist in the vicinity of the Larmor frequency, could be important for applications. 
Therefore, in this paper we present a detailed analytical and numerical analysis of the above 
mentioned problem. 

The paper is organized as follows. In Sec. Ull we describe the model and define the 
lifetime of both the P and Q modes. Here we also derive the boundary conditions and 
transformation properties of the lifetime. In Sec. IHIl we develop an analytical method for 
calculating the frequency dependence of the lifetime of the P mode in the case of high 
anisotropy barrier. Our numerical results obtained by the simulation of the deterministic 
and stochastic Landau-Lifshitz equations are presented in Sec. IIVI Specifically, the features 
of the P and Q modes at zero temperature are studied in Sec. IIV Al and the effects of thermal 
fluctuations are considered in Sec. IIVBI Finally, in Sec. |V]we summarize our findings. 



II. LIFETIMES OF THE PRECESSIONAL MODES 



A. Basic equations of the model 

To study the influence of the rotating magnetic field on the thermal stability of the 
nanoparticle magnetization, we use a minimal model with coherent spin dynamics. Within 
this model, which is applicable to particles whose exchange energy is comparatively large, 
the magnetic state of each particle is completely characterized by the magnetic moment 
m = m(t) of a fixed magnitude m = |m|. Due to its interaction with a heat bath, m(t) is a 
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vector random process which can be described by the stochastic Landau-Lifshitz equation 

-f-m = —7m x (H + n) — — m x (m x H), (2.1) 
at in 

where 7(> 0) is the gyromagnetic ratio, A(> 0) is the damping parameter, the cross denotes 

the vector product, and H = H(t) and n = n(t) are the effective magnetic fields. The first 

field is given by H = — dW/ dm, where W is the magnetic energy of the particle which, in the 

case under consideration, contains only the uniaxial anisotropy energy (l/2)H a m(l—m 2 z /m 2 ) 

and the Zeeman energy — m ■ h(t), i.e., H = H a (m z /m)e z + h(t). Here, the z-axis of a 

Cartesian coordinate system xyz with unit vectors e x , e y and e z is chosen to be parallel to 

the easy axis of magnetization, H a (> 0) is the anisotropy field, m z = m- e z , the dot denotes 

the scalar product, and h(t) is the rotating magnetic field. We assume that h(t) is applied 

perpendicular to the z-axis, so that 

h(i) = h cos(ut)e x + p^isin^^e,,, (2.2) 

where h = \h(t) \ is the field amplitude, u> is the angular rotation frequency, and p = — 1 (for 
clockwise rotation) or +1 (for counterclockwise rotation). It is this field that induces the 
precessional modes of m(t). 

The effective field n(t) accounts for the thermal fluctuations. It is assumed that the 
Cartesian components n a (t) (a = x,y,z) of n(t) are independent Gaussian white noises 
with zero means and correlation functions (n a (ti)n a (t2)) = 2A5(t 2 — ti). Here, the angular 
brackets denote averaging over all realizations of n(t), A = XksT '/'-fin is the noise intensity, 
ks is the Boltzmann constant, T is the absolute temperature, and 5(t) is the Dirac S function. 
In accordance with this definition of n(£), the random process m(t) is Markovian and can 
be described within the Fokker-Planck formalism. 

Because Eq. (12.11) preserves the length of the magnetic moment m(t), it is convenient 
to write the Fokker-Planck equation that corresponds to Eq. (12. ip in spherical coordinates. 
Introducing the polar and azimuthal angles 6(t) and ip(t) of m(t), the forward and backward 
Fokker-Planck equations for the conditional probability density P = P(6,ip,t\6' ,t') (t > 
t') can be written as^ 1 ^ 

d 2p 1 d 2 P d ( „ 2a .„ ,\ „ 



89 2 sin 2 6 dip 2 89 V A 



and 

d 2 P 1 d 2 P ( n , 2a ,..\dP 

2a . „. ,dP 2a dP , 
+ - W ^)- M - = -- w (2.4) 

respectively. Here, t = u r t is the dimensionless time, u r = ^B. a is the Larmor frequency, 
Co = oo/u r is the dimensionless frequency of the rotating field, and a = mH a /2kBT is a 
dimensionless parameter that characterizes the anisotropy barrier height in the units of the 
thermal energy fc#T. Finally, the variables 9 and ib are associated with 9(t) and xb(i) = 
(p(t) — put, respectively, and the functions u(9,ip) and v(9,ip) are expressed through the 
dimensionless magnetic energy 

W = — = - sin 2 9 - h sin 9 cos ip (2.5) 
mH a 2 

(h = h/H a ) as follows: 

= —A sin 9 cos 9 — h sin ip + Xh cos 9 cos ip, 
1 / . „ d . d 



t , (M) = _( sm() __ A _) 



W 



~ ~ sin ib 

= cos9 — hcot9 cosip — Xh -. (2.6) 

sm 9 

It is assumed that the probability density P satisfies the initial condition P\i = ii = 5(9 — 
9')5(ip — ib'). Moreover, if the absorbing boundary conditions are not imposed, then P is 
properly normalized: J 27r dib Jq d9P = 1. 

The system of two stochastic differential equations 

±£{t) = u(9(t),m) + Ya COt9{i) + r?l(f) ' 

libit) = v(9(t),m)-pQ+\^-^- rV2 (t), (2.7) 
at \ a smfc'(t) 

where f]j{t) (j = 1,2) denote independent Gaussian white noises with zero means and 
correlation functions (^(^^(^l)) = 5(t 2 — ti), leads to the same Fokker-Planck equation 
( I2.3p . 33 This means that the above system is stochastically equivalent to the stochastic 
Landau-Lifshitz equation (12. ip . In what follows, we will use Eq. (12. 7p to numerically study 
the thermal stability of the precessional modes of the magnetic moment m(t). 



B. Definition, boundary conditions, and transformation properties of the lifetime 

At zero noise intensity, the rotating magnetic field h(£) can induce stable precessional 
modes of m(t) of two types.— >^ In the first, P mode, the precession angle 0(t) is a constant 
and m(t) in the laboratory frame is a periodic function of time. In the second, Q mode, 
the precession angle varies periodically and m(t) becomes a quasi-periodic function of time 
(because the periods of 0(£) and h(t) are in general not commensurable). Some properties 
of these modes related to the steady-state will be considered in Sec. IIV A[ Here, we use only 
the fact that, depending on the parameters h,u,p, and A, one or two precessional modes 
may exist in steady state. In the latter case, one mode occurs in the up state (a = +1) and 
the other in the down state (a = —1) of m(t). We assume that for a given set of the above 
parameters the magnetic moment m(t) is in the state a if m z (t) tends to am as h slowly 
decreases to zero. It should be noted that the last condition is important because a sharp 
decrease of h can switch m(t) to another state. 

In steady-state, the reference modes are stable and transitions between them are impos- 
sible. However, these transitions can occur under thermal fluctuations. In this case the 
precessional modes become metastable and the magnetic moment remains in a given state 
a for some (dimensionless) random time t a . The average value of this time, i.e., the lifetime 
T a of the metastable state, can be determined using the MFPT method.- 1 ^ The basis of 
this method is the backward Fokker-Planck equation (12. 4p . which should be written for a 
given state a. To this end, we add the index a to all angle variables, replace the conditional 
probability density P by P a = P a (9 (7 ,'ip a ,t\9' a .,ip' a ,t'), and assume that 9 + i,9' +1 E (0,9 ), 
9_ 1 ,9'_ 1 e (-7T - 9 ,n) and i) a ^' a G (0,27r). Here, the angle 9 (n/2 < 9 < n) should 
be chosen so that the time average of the precession angle CT (t) [we recall that Q a {t) for 
the P modes does not depend on t] is relatively close to 9q and to tt — 9 for o = — 1 and 
+1, respectively. To meet these requirements, in our numerical simulations we assume that 
9 = 0.8tt. 

In accordance with the MFPT approach, we consider a circular cone with the cone angle 
7r(l — er)/2 + <r#o as the absorbing boundary for the magnetic moment in the a state, i.e., 
P a \e , (7 =TT(i-a)/2+a9o = 0. In this case, taking into account that P a = P a (9 a ,ipa,u\9' a ,4' / a ,0) 
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(u = t — t') , the lifetime can be denned as T a = J °° duQ a , where 

Qa = # CT / d9 a P a (2.8) 

JO J(n-do)(l-a)/2 

is the probability that the magnetic moment stays in the state a up to a given value of the 
time difference u. Finally, using the relations Q a \ u =o = 1 and Q a \u=oo = 0, one can make 
sure that the lifetime T a = T a (0' a ,il>' a ) is governed by the partial differential equation 

d 2 % 1 d 2 % ( 2a \d% 

= (2.9) 
At the absorbing boundary, the solution of this equation must satisfy the condition 

7^l^=7r(l-CT)/2+<T0 o = 0- (2-10) 

One more important property of the lifetime is that it is a finite function of Q' a and ip' a . 
In order to prove this statement, let us first approximate the stochastic dynamics of the 
magnetic moment by a random walk on the sphere characterized by a dimensionless discrete 
time t = nr, where n = 0, 1, . . . and r is the time step. Then, denoting r an the probability 
that the magnetic moment stays in the state a after the nth step, we can write Q a = 
Tln=i r cm- If the maximal element of the set {r an } equals R a then Q a < R^ T and, as a 
consequence, T a < J °° duR^ T = r/\hiR a \. Finally, taking into account that the condition 
9 < 7r implies that R a < 1, we obtain the desired result: T a < oo. It should also be 
noted that since the maximum angular distance to the absorbing boundary occurs at Q' a = 
7r(l — er)/2, i.e., max 7^- = %\ e , =7r (i_a-)/2' ^ ne condition of finiteness of the lifetime can be 
written in the form 

7^1^=^(1-^/2 < 00 ■ (2- 11 ) 

The above result shows the importance of knowing the solution of Eq. ( 12.91) in a small 
vicinity of the point Q' a = 7r(l — cr)/2. Assuming that Q' a = tt(1 — cr)/2 + > 0), this 

equation at £ CT — > reduces to 

d 2 % 1 d 2 % 10% 2ah, 

+Xs ^wr- 2 i- (2 - 12) 



Its solution can be represented as 

oo 

7; = cln£ 7 + /o + E/KV/)^, (2.13) 
1=1 

where c and fo are constant parameters, and the functions fi = fi(ip' a ) satisfy the ordinary 
differential equations d 2 fi/difj' 2 + /i = and 

yg+Ui- — (-cos^ + Xsm^)— = --5 l2 (2.14) 

with / > 2 and 5i n being the Kronecker delta. According to Eq. (I2.13p . the finiteness 
condition ( 12. lip takes the form c = 0, which in turn is equivalent to lim^^o iadT a /d^ a = 0. 
We note that at c = the derivative &T a /d9' (T in the point Q' a = 7r(l — cr)/2, in general, does 
not vanish: <97^-/96 l ^|^ =7r (i_ (J )/2 = ufi- However, if h — then f\ = (this is so because 
in this case T a does not depend on ■?//.), and the finiteness condition ( 12.1 ip reduces to the 
reflecting boundary condition 

° Tn 0. (2.15) 



d9' a e; =7r (i- CT )/2 

It should also be noted that at h ^ the same reflecting boundary condition holds for the 
average lifetime T a = (l/27r) J 27r dip' a T a , since f 1 = 0. 

Now, using the general equation (12.91) . we can establish the transformation properties of 
its solution T a = T a (9 / a .,ip' a , p) (for clarity, the dependence of T a on p is shown explicitly), 
which is assumed to obey the conditions ( j2 .101) and ( 12. lip . Toward this end, let us introduce 
the change of variables 

£ = ti, = 2*-1/-.. (2-16) 

Taking into account that u{Q' a ,ty'^) = —u(9 f _ a ,ij/_ a ) and u{9' a1 ip' a ) = —u(9'_ a ,ilj'_ a ), one can 
make sure that Eq. ( 12. 9 p in the new variables 9'_ a and ip'_ a becomes 
d 2 H a ^ 1 d 2 U a ( , „, 2a 



89% sin 2 9'_ a d^ ,2 a 



cot^ + y^-.,^)) 



(2.17) 

dH a 2a „ dU a 2a 

where 7i a = T a {ji—9'_ u , 2n—ip / _ a , p). In accordance with the transforms ( 12. 16ft . the conditions 
( 12.1 0p and ( 12. lip for H a take the form 'H (T |e^ (T =7r(i+ CT )/2- CT 0o = and %r|0'_ CT =7r(i+o-)/2 < oo, 
respectively. Therefore, comparing Eqs. ( 12 .9p and ( 12. 17ft and the corresponding absorbing 
and finiteness conditions, one can conclude that H a is equal to T- a (9'_ a ,i/;'_ a ] — p), i.e., 

%(9' a , fa p) = T- CT (7T - 9' a , 2vr - fa -p). (2.18) 



For the average lifetime TV this transformation property reads: T c {§' a 'i P) = T^-a^ — O^; —p). 

In general, the arguments Q' a and ip' a of T a are arbitrary and can be properly chosen to 
best suit the problem. In this paper we are interested in the lifetime of the precessional 
modes reaching the steady state. Therefore, the angles Q' a and ip' a should be associated with 
the solution of deterministic (when a = oo) Landau-Lifshitz equations (12. 7p at some time 
t = t st , i.e., 9' a = 8 a (t at ) and ip' a = Vvfe)- To be sure that the magnetic moment is near the 
steady state, we assume that t st > to + trei, where to is the time of increasing the rotating 
field amplitude from to a given value h, and t re i = 2/ A is the relaxation time. It should 
be emphasized that the time to must be chosen large enough to prevent the dynamical 
switching from the state a to the state —a. In the P mode, the angles 6 a (t s t) and i/j a {tst) 
approach the limiting precession angle CT = lim^^ Q a (t) and the limiting difference of 
phases ^ a = lim^^ ^/ CT (t), respectively. As a consequence, the lifetime of this mode 

r CT = r CT (e CT ,^ ; p) (2.19) 

does not depend on the precise choice of t st . In contrast, in the Q mode the precession angle 
CT (t) is a periodic function of time t with a period Tq and the difference of phases is given 
by ty a (t) = —vt + ^ a (t), where v > and $ CT (t) is also a periodic function with the same 
period Tq (see Sec. IIV Ah . Thus the lifetime of the Q mode, in general, depends on t st : 

% = Ta{e a {t st ), V a {t st ); p). (2.20) 

However, at a ^> 1 this dependence is very weak and can be safely neglected (see Sec. HVBl) . 

In the case of P mode the limiting angles CT and depend also on the parameters 
h, u, X and p. However, for brevity, we keep only the parameter p, i.e., Q a = Q a [p) 
and = ^o-fp]? which together with the state parameter a describes the transformation 
properties of these angles. To find them, we use Eq. ( 12.61) for representing the stationary 
Landau-Lifshitz equations u(Q a [p], ^ a [p\) = and v (Q a [p), ^ CT [p]) = pu as 

cos^ CT (p) = ^(sinQ^p] - pKtanO CT [p]), 
h 

sin^H = -^sinG^p], (2.21) 
h 

where k — + A 2 ). From these equations, it is straightforward to obtain the desired 
result: 

e a [ P ] = vr - e_ CT [-p], * a \p] = 2vr - *_ CT [-p]. (2.22) 
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Because the transformation properties (12.221) are similar to those in Eq. (I2.16p . from 
Eq. (I2.18P one gets the transformation property of the lifetime of the P mode 

%(Qa[p\^a[p\; p) = T-^tt - CT [p], 2tt - * a [p\; -p). (2.23) 

It shows that the lifetimes characterized by the pairs {a, p} and {— er, — p} are the same. 

III. ANALYTICAL SOLUTION OF EQ. (^9j) 
A. Three-mode approximation 

The analytical determination of the lifetimes of the precessional modes implies the solu- 
tion of Eq. (12. 9 p with the absorbing and finiteness conditions (I2.10p and (12. lip . Since the 
lifetime T a is a periodic function of ip' a with the period 27r, it can be expressed as the Fourier 
series 

oo 

T = £ % n (9' a )e m ^. (3.1) 

n=— oo 

To guarantee the reality of T a , we assume that the coefficients T an = Tan{9' a ) of the series 
satisfy the condition T a - n = T* n (the asterisk denotes complex conjugation). Substituting 
this series into Eq. f 1 2 . 9 1) and introducing the differential operators 



d 2 , „, . s d ii 



2 



2 , 



L ^=d^ + ^t9' a -asin29' a )—- ^ ^ 

flan. n . „. 
+ i—{cos9 a - pu), 

N n =(\co S 9 , a -i)^ + -^--incot9 , a , (3.2) 
d9 a sm v' a 

one obtains for T an an infinite set of coupled equations 

L n Tan + —^{.Nn+iTan+l + NZ n+1 Tan-l) = — — 5 n Q. (3-3) 

Like T a , the coefficients Tan must satisfy both the absorbing boundary condition 
Tan\e' a =iv{i-a)/2+ae = and the finiteness condition T an 1^=^(1-^/2 < oo. 

In the three- mode approximation, when T an — for all \n\ > 2, an infinite set of equations 
(13. 3 p reduces to a set of three equations for Ta®, T a i and T* v Since T a is real, it is convenient 
to consider instead of T a \ and T*i the real and imaginary parts of T a i, i.e., T<X = KeTai and 
Tai = Im7^i. Then, in this approximation, the lifetime (13. ip is given by 

Ta = Tao + 2T+ cos ^ - 2T a \ sin ^ (3.4) 
10 



and, denoting the real and imaginary parts of an operator O as + and 6 , respectively, 
from Eqs. (13 .3p and (13. 2p one obtains a set of three coupled equations for T^o, TX an d T^i- 

LoXo + y (1 + hNt%\ - hN^T-,) = 0, (3.5) 
LfT+ - L- X r- X = —NfToo, (3.6) 
Ltr~ x + L- X T+ = yiV -T CT0 . (3.7) 

Using Eq. (I3.5p . the function 7^o(^) can be expressed through TXi^'a) an d T&iiQ'o)- 
Indeed, considering 

FM) = hNtXX - hNrXi (3.8) 
as a given function of 6' a , one can write a formal solution of Eq. (I3.5p .^ From this solution, 
by satisfying the absorbing and reflecting boundary conditions 7^o|^=7r(i-o-)/2+o-e = an d 
dTa-o/d9' a \e' cT = 7T (i-a)/2 = (since T a o does not depend on 9' a , the last condition is equivalent 
to the finiteness condition T a o\e' a =n(i-a)/2 < oo, see Sec. Ill B[) . we obtain 

2a /^cos6»; e -ax 2 r i 
T a0 = — dx- / dy[l + F CT (arccos ay)]e ay . (3.9) 

A Jcos 8o 1 — X J x 

Substituting this result into Eqs. (13 .6p and (13. 7p . one can readily get the coupled integro- 
differential equations for and T~\ which, however, are too complicated to be solved in 
the general case. Moreover, these equations are not closed because, in accordance with 
Eqs. (13. 8p . (13. 6 p and (13. 7p . the function F CT (arccos<7?/) depends on 7^-o- Fortunately, an 
important case of high anisotropy barrier can be studied in detail. 



B. High anisotropy barrier 

In the case of high anisotropy barrier (or low temperatures), when the condition a = 
mH a /2kBT ^> 1 holds, the main contribution to the first integral in Eq. (13. 9p comes either 
from a small vicinity of the point x = (if Q' a is not too close to tt/2 and acos^ > 0) 
or from the vicinity of the point x = a cos Q' a (if <rcos6^. < 0). It should be noted that 
the latter situation can be realized only in the case of Q mode with max9 + i(t) > 7r/2 or 
min B_i(t) < 7r/2. Here, we restrict our theoretical analysis by considering small amplitudes 
of the rotating field that do not exceed the threshold amplitude of the Q mode. Such rotating 
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field can induce only the P mode and, as a consequence, the former situation is always 
realized (see also Sec. IIV Ap . Therefore, taking into account that in this case the integrals 



/cos 0o e " ^ xe ^ an d Jo dye ay2 at a > 1 can be approximated by /f^ dxe ax2 = \J^/a and 
e a /2a, respectively, Eq. f!3 . QH reduces to 

% = ^ x f^ (l + 2ae- a f 1 dyFJaiccosay)e ay2 \ . (3.10) 
A V a V Jo J 

In order to evaluate the integral in Eq. (I3.10p . we note that small vicinities of the limits 
of the integration determine the asymptotic behavior of this integral at a — > oo. More 
precisely, the lower limit is responsible for the high-frequency behavior of this integral and 
the upper one for its behavior in the resonant case. Hereafter, we call the rotating field 
resonant if cD ~ 1 and the direction of its rotation coincides with the direction of the natural 
precession of the magnetic moment, i.e., if ap = +1. To study these cases in a unified way, 
it is convenient to write the integral in Eq. (I3.10p as a sum of two terms, F a (&rccos a)e a /2a 
and Fo-(arccos 0), which correspond to the resonant and high frequencies, respectively. Thus, 
taking into account that arccosa = 7r(l — a)/2 and arccosO = ir/2, Eq. (I3.10p yields 



%0 ~ Tvf 



(3.11) 



It is important to stress at this point that Eq. (13. lip is not an exact asymptotic formula for 
To-o because the terms ^(^(l — er)/2) and 2ae~ a F a (7i/2) correspond to different frequencies. 



1. Vicinity of the point 6' a = n{l — o~)/2 

To calculate F CT (7r(l — cr)/2), we need to solve Eqs. (13.61) and (13.71) in the vicinity of the 
point 9' a = 7r(l — ct)/2. Assuming that 9' a = 7r(l — cr)/2 + arj a (0 < T] a <C 1) and ah 1, 
these equations can be rewritten as 

LtT^-^a-p^T-^-ah^-, 



where 

d 2 ( 1 \ d 1 , 
L+ = — +(--2a Va )- -. 3.13 
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Since the approximate formula (13.1 ip does not depend on 9' a , it cannot be used to determine 
the derivative dTao/dr/^. Therefore we use an exact result: 

= -^OrL^ 1 + ^—y*"> • (3.14) 
following from Eq. ( I3.9p . which at r\ a -C 1 gives 

dTao a 



l + /v(|(l-a; 



Vv (3.15) 



It is not difficult to verify that an exact solution of Eq. (I3.12p that vanishes as h — > has 
the form 

(1 — A 2 — opuS)h dTao 



'crl 



2[(1 -<rpu;) 2 + A 2 ] dr/ ff 



_ A (2 - <rpu;)/i rf^o 
^"^[(l-apc^ + A 2 ] ^' ^ ibj 

where dT <T o/dr] a is given by Eq. ( I3.15p . According to this equation, the solution ( I3.16P 
contains an unknown parameter F a (ir(l — a)/2), which can be determined from the fitting 
condition F a {-n{\ — er)/2) = (hN^T^ — hNiT~i)\ v „=o [see Eq. (13. 8p ], Taking into account 
that Ni = Xd/drjcr + X/rj^ and iV-f = —odjdri a — cr/rj a , this condition reduces to 

Finally, substituting Eq. (I3.15P into Eq. ( I3.17P and imposing the commonly used condition 
A 2 1, one gets 

/7T. A ah 2 

F a -(1-0-) ] = ^. 3.18 

V2 V 7 (1 - apu) 2 + A 2 + ah 2 V ; 

2. Vicinity of the point Q' a = tt/2 

To find F CT (7r/2), we assume that 6' a = tt/2 + £ a (|£ CT | <C 1) and, as before, ah <C 1. Since 
in this case L{ = — 2a(£ CT + puj)/X, Nq = —X^ a d/d^ a and Nq = —d/d^ a , Eqs. (I3.6P and 
(13. 7p take the form 

L tr- 1 - 2 ^ + ^)TA = -^, (3.19) 
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where 



d 2 n , d 



LJ = — - 2(/^~ 1 

and, as it follows from Eq. (I3.14p . the derivative dTa /d£ a at £ CT — > is given by 



d% 



crO 



di„ 



-a- 



X 



71 



l + FJ-(l-a)\ +2ae- a FJ- 
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(3.20) 



(3.21) 



In general, the solution of Eq. ( 13.19)) can be represented in the form of the Taylor se- 
ries: 7^1 = Y^=o c tnia- However, since the main quantity of our interest is F a {-n/2) = 
(hNi TX ~ hNiT^l^^ = A/ic^ + hc~! (here, = A and iVf = —d/d^ a ), we can restrict 
ourselves to the linear approximation, i.e., 7^ = c^ + c^ x ^ a . Assuming that Cj ^> A/2a, 
after straightforward calculations one obtains 



Xh dT a o 
4ad) 2 d^ a ' 



h dT a0 



and 



A(l + Co 2 )~h dT a0 



2oj d^ a 

h d%o 



(3.22) 



(3.23) 



al H 2Co{X 2 +u 2 ) d£ a ' al " 2{X 2 + Co 2 )d^ 
Therefore, the fitting condition F a {7i/2) = Xhc^ + hc~ x can be reduced to the following one: 

A(l + A 2 + 2u 2 )h 2 d% 



P- 



2Cj(X 2 + Cj 



d^ 



which together with Eq. (I3.2ip at A 2 1 yields 

e a [l + F a (7t(l-a)/2)]h 2 



-op- 



2[g(Q) + apah 2 } 

where g(u) = Cj(X 2 + w 2 )/(l + 2u 2 ) and F ff (7r(l - a)/2) is defined in Eq. fl3TT8|) . 



(3.24) 



(3.25) 



5. Frequency dependence of the lifetime 

Now we are in a position to determine the lifetime ( 12.1 9p of the P mode. Using Eq. (13.4)) 
and the fact that T a o does not depend on 0' a and ip' a , the desired lifetime can be expressed 
as 

% = % + 2{%\ co^m - r; x sin =9aM . (3.26) 

If the field amplitude ft, is small enough, then cos 0c- [p] in Eq. (12.211) can be replaced by a. 
One can easily check that in this case 

***M = -r 1 + f_ rfB , (3^27) 

14 



and the precession angle CT [p] is given by 



7T 



Q a [p] = (l- a ) + ah A 



l + X2 (3.28) 



2 V ' \ (1-opCj) 2 + A 2 ' 

The last result shows that the functions Tji in Eq. ( I3.26|) must be taken from Eq. (13 . 1 6[) 
with r] a = crO CT [p] + 7r(l — o)/2. But, according to Eq. ( 13. Ill) , the term 7^o contains an 
additional factor e a , and so the second term in Eq. (13.261) which describes the dependence 
of T a on ^ a [p] can be safely neglected at a ^> 1. 

Thus, using Eqs. (13. lip . A3. 18j) . and (I3.25p . for the lifetime of the P mode we obtain 

% = j^R ap (uj)S ap (uj), (3.29) 

where 

(« (Xi lb , . 

Cj) = 1 =- 3.30 

' (1 - opCj) 2 + A 2 + ah 2 K ' 

and 

S„ P W = 1 - af>ak2 ~ ■ (3.31) 
g(u) + a pah 2 

If ap = +1 then Eq. (13. 29j) correctly describes the frequency dependence of the lifetime 

in the vicinity of the point Co = 1 and at Cj ^> 1. Since ah 2 <C 1, A 2 1 and g(l) ~ 



1/3, in the former case we obtain T a = (e a /A)y it / aR + i(Co) . To put it differently, the 
rotating magnetic field whose direction of rotation coincides with the direction of the natural 
precession decreases the lifetime in a resonant manner, i.e., a resonant suppression of the 
thermal stability of the P mode occurs. Taking into account that 1 — R + i(Cj) oc Cj" 2 and 



1 — S + i(u) oc Cj 1 , in the latter case Eq. (I3.29j) yields T a = (e a / n/aS + i(u). In contrast, 
at op = —1 Eq. (13.291) describes only the high-frequency behavior of the lifetime: T a = 



e a /\)\Jn/aS-i(uj). Combining the last two results, we obtain the expression 

_ e a fW ( 2ah 2 \ 

^-ti: 1 -^ ( 3 - 32 ) 



A V a V Cj 

(Cj ^> 1), which shows that while at up = +1 the rotating field suppresses the lifetime of the 
P mode, the rotating field with op = — 1 enhances it. 



4- Lifetime at zero frequency 



For the evaluation of the lifetime at u = it is convenient to use the one-dimensional 
approximation, which consists in replacing ip by in the magnetic energy (12. 5p . In this 
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case u(9' a ,ip' a ) = — (A/2) sin 26^. + Xhcos9' a , the lifetime T a does not depend on ip' a , and the 
partial differential equation (12. 9 p reduces to the ordinary one: 

^ + (cot a - a sin 2ff a + 2ah cos O^f = ~ ■ (3-33) 

Its exact solution satisfying the conditions ( I2.10p and (12. 15ft [we note that for this equation 



the fmiteness condition (12. lip is equivalent to the reflecting boundary condition (I2.15P 
be written in the form 



can 



% = dx\ / dye^\ (3.34) 



A J cos 6*o 1 — X z 

where f(x) = x 2 +2h\/l — x 2 (\x\ < 1) is the symmetric function with min/(x) = /(0) = 2h, 
maxf(x) = f(y/l- h 2 ) = l + h 2 , and f(l) = 1. 

As before, we are interested in the behavior of T a at a ^> 1. In this case a small vicinity 
of the point x = gives the main contribution to the first integral in Eq. (I3.34p . This 
contribution depends not only on a but also on the parameter ah. In particular, if ah <C 1 
then, putting x = everywhere except for e~ a ^ x \ representing e~ a ^ x ^ in the vicinity of the 
point x = as e - 2ah - a ( l - h ) x2 ~ (i _ 2ah)e~ ax2 and extending the limits of integration to 
infinity, Eq. d3~34l) yields 

% = yVtto (1 - 2ah) f 1 dye af{y) . (3.35) 



A ' Jo 

Then, taking into account that 



f 1 dye af{y) « e a [°° dze- 2az (l + 2 3/2 ah^) 
Jo Jo 



, ,1 + Vmh) (3.36) 



(a 3> 1, ah «C 1) and neglecting terms proportional to -s/Trah, from Eqs. (I3.35P and ( I3.36P 
one finds 

%=^J^{l-2ah). (3.37) 
A V a 

Comparing this result with Eq. ( I3.32p . we conclude that 7^-|^=o < T a \u=oc, and so at ap = 
— 1 the frequency dependence of the lifetime has a maximum exceeding the limiting value 
T a \u =OQ = {e a /\)^fnja. 

It should be noted that the lifetime 7^|^ = o strongly decreases with increasing h. For 
example, if h ~ l/\/a then, instead of Eq. ( I3.35p . we obtain 



% = \ x i^re- 2ah I'dye^v). ( 3 .38) 
A y 1 - h Jo 
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To evaluate the integral in Eq. (13. 38 p . it is convenient to divide the interval of integration 



(0, 1) into two parts, (0, y 1 — h 2 ) and (y 1 — h 2 , 1), and apply the Laplace method. 36 By 
this way, the corresponding integrals can easily be evaluated yielding 



Jo \ a Jo 



^e^~ h %rf(^) (3.39) 



and 



2 V a V h 



J\/l-h 2 y/a Jo 

where erf(z) = (2/y/Tr) Jq dxe~ x2 is the error function. Finally, using Eqs. ( I3.38[) ( ]3.40p and 
the approximate formula erf (y/a/h) ~ 1 (y/a/h ~ a 1), we find the following expression 
for the lifetime at Q = 0: 



-z 2 



AVI - h 



2 J 



(3.41) 



which is valid if a 3> 1 and ft, is of the order of l/\fa. It is not difficult to see that for 
these conditions the strong inequality T a \cj=o 7^-1^=00 holds. Taking into account also 
Eq. (13.321) . we can conclude that T a as a function of u) at op = — 1 has a local maximum 
(see Sec. UVBl) . 



IV. NUMERICAL RESULTS 



A. Precessional modes of the magnetic moment 

The analytical results suggest that the numerical analysis of the lifetimes of the preces- 
sional modes should start with the study of these modes without thermal fluctuations. More 
precisely, it is necessary (i) to determine the conditions when for a given rotating field one 
precessional mode exists in the up state of the magnetic moment and one in the down state, 
and (ii) to study the steady-state properties of these modes. To solve these problems, we use 
the Landau-Lifshitz equation written in the form of Eq. (12.71) with a = oo. In general, the 

17 



finite state of the magnetic moment (i.e., state at t — > oo) depends not only on the parame- 
ters h, u and p characterizing the rotating field, but also on how this field is switched on. In 
particular, a sharp switching of the rotating field induces dynamical effects that may result 
in a change of the initial state ex.— Although these effects can be important for applications, 
they are out of our scope here. Therefore, to minimize the role of dynamical effects, we 
assume that the switch-on of the rotating field is slow enough. 

In order to determine the character of the precessional modes for a given rotating field, 
the following numerical scheme is used. First, the field amplitude is discretized as h = nAh, 
where n = 1, 2, . . . and Ah is the amplitude increment. Then, putting n = 1 and using the 
initial conditions #(0) = 7r(l — cr)/2 + alCT 4 and ^(0) = 0, the fourth order Runge-Kutta 
method with a time step At = 1CT 3 is applied to solve Eq. (12. 7p at a = oo on the time 
interval [0, i m ] (t m 3> t Te \ — 2/ A). It is assumed that at i = t m the field amplitude jumps 
from Ah to 2 Ah, i.e., n becomes equal to two, and Eq. (12. 7p is solved again on the interval 
[0,t m ]. But now the initial conditions are the solutions obtained for t = t m at the previous 
stage: 0(0)\ n=2 = #(t m )|n=i and ip(0)\ n= 2 = V^nOlr^i- Continuing this procedure, Eq. (I2.7P 
can be solved for an arbitrary n. It is worth to note that since t m 3> t re i the solutions of this 
equation at i ~ t m are expected to be quite close to the steady-state solutions Q a (t) and 

Using the above procedure with Ah = 10 -2 , t m = 10 3 and A = 0.15 [this value of A, 
which belongs to the interval (0.01, 0.22) of typical values of the damping parameter in the 
case of Co samples, 37 is used in all our numerical calculations], we determined the character 
of the precessional modes for a wide range of parameters characterizing the rotating field. 
It is established that if ap — — 1 then only the P mode is realized for all h and u. In 
contrast, the precessional modes at ap = +1 exhibit a much more complex behavior. The 
results related to the character of these modes are summarized in the diagram shown in 
Fig. [TJ We note that the difference between two P modes, which exist in the regions P + i 
and P+x, is that the precession angle Q a as a function of h is discontinuous at the boundary 
between them. It should also be emphasized that the transitions between the modes with 
ap = +1, which occur under changing the field amplitude h, are reversible. For clearness 
of presentation, this fact is illustrated by the horizontal bidirectional arrows. In contrast, 
the transitions to the P mode with ap = — 1 are irreversible (they are depicted by the 
horizontal unidirectional arrows). Moreover, using Eq. (I2.2ip and the stability criterion for 
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FIG. 1: (Color online) Diagram of the precessional modes for ap = +1. The regions in the h-Co 
plane where different P modes exist at ap = +1 are denoted as P+i (white) and Ptj (light-green). 
The Q mode is realized in the white shaded region. In the region denoted as P_i (blue) the 
stable precessional modes with ap = +1 do not exist. Here, only the P mode with ap = —1 is 
realized. The vertical dotted lines (a), (b), (c) and (d) correspond to h = 0.05, 0.1, 0.18 and 0.25, 
respectively. 

the P mode,— we independently confirmed the correctness of this diagram by calculating 
the lines that separate the regions with ap = +1. 

Since the rotating field is switched on during the time interval of duration to — (n — 1)10 3 , 
the above procedure is time consuming. However, in view of the importance of the diagram 
of the precessional modes for the problem of lifetimes, its use for the precise determination 
of this diagram is quite acceptable. At the same time, the application of this method 
to the study of the steady-state properties of a given mode, whose character is already 
known from the diagram, is clearly redundant. Therefore, to reduce the computational 
time, next we use a modified numerical procedure with t m = At and Ah = hAt/50 leading 
to to = Ath/Ah = 50. 

Figure [2] shows the frequency dependence of the precession angle G+i for different pre- 
cessional modes. The dashed (black) and solid (blue) lines correspond to the P modes with 
ap = — 1 and ap = +1, respectively. In accordance with Fig. [TJ the Q mode occurs at 
cl>2 < uj < 0J3 (we recall that this mode can exist only if ap = +1). For this mode, the 
time dependence of the precession angle 0+i(t) and the time dependence of the difference of 
phases ^+i(t) = —vt + $+i(t) are illustrated in Figs. [3]and|U We note that for a given set 
of parameters max0 + i(t) > tt/2, i.e., the precession angle Q a (t) in the case of Q mode can 
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FIG. 2: (Color online) Frequency dependence of the precession angle G+i for different modes that 
exist at h = 0.25. The frequencies u\ = 0.49, 0J2 = 0.70 and 003 = 0.89 are the coordinates of the 
points in which the vertical dotted line (d), see Fig. [IJ crosses the boundaries of the diagram (cl>i, 
C02 , and 0)3 depend on h). The green line (with squares) and brown line (with triangles) show the 
frequency dependence of maxG + i(i) and minG_|_i(t), respectively, in the case of Q mode. 



1.8 




FIG. 3: (Color online) Time dependence of the precession angle ® + \(t) in the case of Q mode. 
The parameters of the rotating field are as follows: p = +1, Co = 0.725, and h = 0.25. 

cross the anisotropy barrier. Moreover, since sin0 CT (t) and sin\l/ -(t) have the same period 
Tq, the parameter v and the period Tq (its frequency dependence illustrates Figs. [5]) are 
connected by the condition v = 2/c7t/Tq, where k is a nonnegative integer that depends on h 
and u. In particular, if h = 0.25 then k = 1 at 0)2 < & < &' 2 , where Q' 2 = 0.78 is the solution 
of the equation min0 + i(t)| (i=(I ,^ = (see Fig. [2]), and k = at u' 2 < Co < 0)3. It should also 
be stressed that in the last case the function $+i(t) shows a cosine-like behavior, in contrast 
to that shown in Fig. HI 
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FIG. 4: (Color online) Time dependence of the difference of phases = —vt + <I>+i(i) in the 

case of Q mode. Insert: time dependence of the function $ + i(i). The parameters of the rotating 
field are the same as in Fig. [3] and v = 0.38. 

To 
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FIG. 5: (Color online) Frequency dependence of the period Tq of the precession angle + i(f) in 
the case of Q mode. The rotating field parameters are taken as p = +1 and h = 0.25. 

B. Simulated lifetimes and their properties 

As it was mentioned earlier, the thermal fluctuations can cause the transitions between 
different precessional modes induced by a given rotating field. According to the diagram in 
Fig. [1} these modes, one in the up state of the magnetic moment and the other in the down 
state, exist only if the amplitude and frequency of the rotating field belong to the region P + i, 
P+i or Q. In this case, the lifetime of a given mode can be calculated from the numerical 
solution of the stochastic equations (\2.7\i . Our testing calculations showed that the solution 
of these equations by the Euler method gives practically the same lifetime obtained by the 
fourth order Runge-Kutta method. But in the first case the calculation time is almost four 
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FIG. 6: (Color online) Schematic time dependence of the polar angle 9{t) in the regions P+i and 
P +1 shown in Fig. HJ The change of the magnetic moment state a from +1 to —1 occurs at 
t = to + i re i + i+i, when 6(t) reaches the angle 9q = 0.8n (the horizontal dashed line) for the first 
time. For a given trajectory 6(t), the lifetime of the P mode in the state a = +1 is equal to t + \. 
Running N 3> 1 trajectories, the mean lifetime can be evaluated as T+i = (l/N) J2?=i t+i- 

times less. Therefore, because the procedure of determining the frequency dependence of 
the lifetime is extremely time-consuming, we used the Euler method. The time step At is 
chosen to be 10~ 3 and the initial conditions are given by 6 a (0) = 7r(l — er)/2 + al0~ 2 and 
"0o-(O) = 0. To prevent the appearance of singularities in Eq. (12.71) at 6 a = 7r(l — er)/2, we 
assume that the point Q a = 7r(l — cr)/2 + erlO -3 acts on the process 9 a (t) as a reflecting 
screen. Since our interest here is the lifetimes of the precessional modes reaching the steady 
state, the thermal fluctuations are switched on Out t — (see Sec. Ill Bp with t st = to + ^rei ( see 
Fig. [6]). In the case of P modes we chose t st = 10 2 , while for the Q mode t st G [10 2 , 10 2 + Tq]. 
In the latter case, the thermal fluctuations can be switched on at a certain instant of time, 
e.g., when the precession angle Q a {t) reaches maximum or minimum. Finally, in all our 
numerical simulations the parameter a = mH a /2kBT is chosen to be ten. 

In Fig. [7J we show the frequency dependencies of the lifetime 7+i for the rotating field 
amplitudes indicated in Fig. [TJ Each point of these curves is determined by running N = 10 4 
trajectories of the polar angle 9{t) (see the caption to Fig.E]). If h < 0.19 then the dependence 
of 7+i on Co is continuous and exhibits a resonant minimum at Cj = uj res . The plot of the 
resonant frequency a) res versus the field amplitude h is shown in Fig. [HJ In contrast, if 
h > 0.19 then 7+i is discontinuous: at Co G (0)1,0)2) the function 7+i does not exist. This 
result is a consequence of the fact that in the region P_i (see Fig. [[]) there are no stable 
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FIG. 7: (Color online) Frequency dependencies of the lifetime of the precessional modes induced 
by the rotating field with p = +1 in the up state (a = +1) of the magnetic moment. 
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FIG. 8: (Color online) Dependence of the resonant frequency of the lifetime 7+i on the rotating 
field amplitude. 

precessional modes with op = +1. 

One more important feature of the lifetime 7+i is that it is practically not sensitive to 
changing the character of the precessional modes. In particular, for h = 0.1 the frequency 
dependence of 7+i (see the red line with circles in Fig. [7j) is continuous at the point u) = 0.74 
(i.e., point separating the regions P + i and at h = 0.1), while the precession angle 
@ + i is discontinuous. This insensibility of the lifetime to changing the precessional modes 
with changing the field frequency Q is especially surprising when the Q mode appears. For 
example, even at h — 0.25, when the precession angle can cross the anisotropy barrier (see 
Fig. [3]), the character of the frequency dependence of 7+i (see the green line with stars 
in Fig. [7]) is changed at w = w 3 so small that it is not visible on this scale. Moreover, 
the lifetime of the Q mode does almost not depend on t st G [10 2 , 10 2 + Tq]. This result is 
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FIG. 9: (Color online) Probability density functions of the lifetime of the Q mode for two values 
of t s t- The blue solid and red dashed lines correspond to such t st that 0+i(i s t) = maxO+i(£) and 
O+i (i s t) = min0 + i(t), respectively. Insert: the same density functions in a larger time scale. The 
parameters of the rotating field are as follows: p = +1, ui = 0.75, and h = 0.25. 

counterintuitive because different t st may correspond to very different values of the precession 
angle + i(i st ). Therefore, in order to get more insight into the problem, we calculated the 
probability density function D(t + i) of the lifetime of the Q mode for two special values of 
t s t which correspond to minG + i(t) and maxO+i(i), respectively. As seen from Fig. [9j these 
density functions are somewhat different from each other only if t +1 < t ie \. In this region 
the probability density D(t + x) depends on t st and exhibits local minima and maxima which 
come from a complex behavior of the Q mode. In contrast, at t+i ^> t Ie \ the memory about 
the chosen value of t st and periodicity of 0+i(t) and $+i(t) is lost. As a consequence, in this 
region the difference between density functions vanishes and the local minima and maxima 
disappear. Such behavior of D(t + i) on t st confirms that the lifetime 7+i = Jo°° t+iD(t+i) 
practically does not depend on t st . 

Finally, the influence of the direction of field rotation on the frequency dependence of the 
lifetime of the precessional modes is illustrated in Fig. [TUJ In accordance with our analytical 
results, the rotating field with p = +1 and p = — 1 influences the lifetime 7+i in a different 
way. Specifically, 7+i as a function of Cj at p = +1 displays a deep minimum, while at 
p = — 1 it shows a pronounced maximum. This difference in the behavior of the lifetime 
results from that the rotating fields with different p act on the magnetic moment in a given 
state a quite differently. From a physical point of view, the reason is that the magnetic 
moment has a definite direction of the natural precession. We note also that the numerical 
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FIG. 10: (Color online) Frequency dependencies of the lifetime of the precessional modes in the up 
state (cr = +1) of the magnetic moment. It is assumed that h = 0.25 for both clockwise (p = — 1) 
and counterclockwise (p = +1) rotation of the magnetic field. 

data 7+i|tD_Kx) = 9.1 x 10 4 and 7+i|^o = 1-97 x 10 3 are in a good agreement with the 
analytical results 8.2 x 10 4 and 2.38 x 10 3 obtained from the asymptotic formulas (I3.32p and 
(I3.4ip . respectively. Some difference between them can be caused by that the asymptotic 
formulas, which were obtained at a — > oo, are applied to a = 10. 

V. CONCLUSIONS 

We have studied in detail the thermal stability of the precessional modes of the nanopar- 
ticle magnetic moment induced by the rotating magnetic field whose plane of rotation is 
perpendicular to the easy axis of the nanoparticle. If the direction of field rotation and the 
direction of the natural precession of the magnetic moment are opposite, i.e., if the condition 
ap = —1 holds, then only periodic (P) stable mode is induced by this field. In contrast, 
if the above mentioned directions coincide, i.e., if ap = +1, then the magnetic moment 
exhibits a much more complicated behavior. The numerical solution of the deterministic 
Landau-Lifshitz equation in the long-time limit has shown that, depending on the rotating 
field amplitude and frequency, in this case the magnetic moment can be in one of two P 
modes, in the quasi-periodic (Q) mode, or even be unstable. These results obtained in the 
absence of thermal fluctuations have been collected in the diagram shown in Fig. 1. 

If the amplitude and frequency of the rotating field are chosen so that a stable precessional 
mode exists in both up (cr = +1) and down (cr = —1) states of the magnetic moment, 
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then the thermal fluctuations can cause transitions between these modes. One of the most 
important parameters characterizing these transitions is the lifetime of a given mode. Since 
it can be naturally associated with the mean first-passage time for the magnetic moment, we 
have used the Fokker-Planck formalism to define this quantity and calculate its properties. 
In particular, we have determined the boundary conditions and transformation properties of 
the lifetime and have developed an analytical method for finding its frequency dependence 
in the case of high anisotropy barrier and small amplitudes of the rotating field. Using 
this method, it has been shown that the rotating field (a) slightly decreases (if ap = +1) 
or increases (if ap — —1) the lifetime of the P mode at large frequencies and (b) strongly 
decreases it (if ap = +1) in the vicinity of the Larmor frequency. We have also established 
that at zero frequency the lifetime is always less than in the limit of large frequencies. 

These analytical findings for the lifetime of the P mode have been confirmed by our 
numerical simulations of the stochastic Landau-Lifshitz equation. Moreover, the numerical 
simulations of this equation for not too small amplitudes of the rotating field permitted us 
to solve the problem of the lifetime of the Q mode. Since in this case the precession angle 
is a periodic function of time which can cross the anisotropy barrier, the solution of this 
problem is of particular interest. It has turned out that, although the precessional angle 
depends on time, the lifetime of the Q mode practically does not depend on this time, i.e., 
on the precession angle. We have also verified this result by calculating the lifetime from the 
first-passage time distributions that correspond to different values of the precession angle. 
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